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Abstract 



A complete classification of the computational complexity of the fixed-point exis- 
tence problem for boolean dynamical systems, i.e., finite discrete dynamical systems 
over the domain {0, 1}, is presented. For function classes T and graph classes G, an 
(jF, t/)-system is a boolean dynamical system such that all local transition functions lie 
. in and the underlying graph lies in Q. Let be a class of boolean functions which 

' is closed under composition and let be a class of graphs which is closed under taking 

O . minors. The following dichotomy theorems are shown: (1) If J- contains the self-dual 

functions and G contains the planar graphs, then the fixed-point existence problem 
. for {!F, C/)-systems with local transition function given by truth-tables is NP-complete; 

^ ' otherwise, it is decidable in polynomial time. (2) If J-^ contains the self-dual functions 

, and G contains the graphs having vertex covers of size one, then the fixed-point exis- 

tence problem for [J-', ty)-systems with local transition function given by formulas or 
circuits is NP-complete; otherwise, it is decidable in polynomial time. 



^ ■ 1 Introduction 

O 

Background on complex systems. A complex system, in a mathematical sense, can be 
viewed as a collection of highly interdependent variables. A discrete dynamical system 
^ ' is a complex system where variables update their values in discrete time. Though the 

■ interdependencies among the variables might have quite simple descriptions on a local 

level, the overall global behavior of the systems can be as complicated as unpredictable or 
undecidable (see, e.g., [El [20]). This phenomenon has been widely studied in the theory 
of cellular automata [3T1|33] and its applications (see, e.g., [Ml [21]). 

Finite discrete dynamical systems are characterized by finite sets of variables which can 
take values from a finite domain. In essence (see, e.g., |101 [9]). a finite discrete dynamical 
system (over a finite domain) consists of (1) a finite undirected graph, where vertices corre- 
spond to variables and edges correspond to an interdependence between the two connected 
variables, (2) a set of local transition functions, one for each vertex, that map values of 
variables depending on the values of all connected variables to new variable values, and 
(3) an update schedule that governs which variables are allowed to update their values 
in a certain time step. A formal definition is given in Sect. [2l Due to their structural 
simplicity and modelling flexibility, finite discrete dynamical systems are suitable for ana- 
lyzing the behavior of real-world complex systems. In fact, the conception is motivated by 
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analysis and simulation issues in traffic flow (see, e.g, [21 E]) and inter-domain routing [18j . 
It also has applications to two-species diffusion-reaction systems such as synchronous or 
asynchronous versions of the nearest neighbor coalescence reaction A + A ^ A on a lattice 
in the immobile reactant case [1]. 

A central problem in the study of discrete dynamical systems is the classification of 
systems according to the predictability of their behavior. In the finite setting, a certain be- 
havioral pattern is considered predictable if it can be decided in polynomial time whether a 
given system will show the behavioral pattern |12] . In a rather strong sense, predictability 
and tractability are identified. It is not surprising that the reachability of patterns is, in 
general, an intractable problem, i.e., at least NP-hard (see, e.g., [171 ETJ Ej). However, 
some restricted subclasses of finite discrete dynamical systems, i.e., systems given by re- 
stricted sets of local transition functions and network topologies, are known to possess 
easy-to-predict patterns (see, e.g., |5l [6l |7j and the discussion of related work below). For 
the purpose of analyzing and simulating real-world systems by finite discrete dynamical 
systems, it is highly desirable to have sharp boundaries between tractable and intractable 
cases. 

A fundamental behavioral pattern for discrete dynamical systems are fixed points (ho- 
mogeneous states, equilibria). A value assignment to the variables of a system is a fixed 
point if the values assigned to the variables are left unchanged if the system updates the 
values. A series of recent papers has been devoted to identification of finite systems with 
tractable/intractable fixed-point analyses [HJ [301 EHl EH] . But although it is an old question 
how common intractability results are for discrete dynamical systems (see \32\ Problem 
19]), a precise characterization of the islands of tractability even for the fixed-point exis- 
tence problem in the simplest case of the boolean domain {0, 1} has remained an open 
problem (see ^). 

Contributions of the paper. We contribute to the problem of classifying boolean (discrete) 
dynamical systems, i.e., finite discrete dynamical systems over the domain {0,1}, with 
regard to the computational complexity of the fixed-point existence problem, i.e., decide 
whether a given system has a fixed point, in two ways. 

A first contribution is the proposal of a general analysis framework for systems. We 
say that a boolean dynamical system is an (.F, t/)-system if its local transition functions 
belong to the class of boolean functions and the underlying graph belongs to the graph 
class G- We propose to consider two well-studied frameworks for functions and graphs (a 
formal introduction is given in Sect. [3|): 

• Post classes [23|: We assume that the function class J- is closed under composition 
(and some more reasonable operations). Equivalently, is a class of boolean func- 
tions that can be built by arbitrary circuits over gates from some finite logical basis 
[23]. Examples are the monotone functions, the linear functions, and the self-dual 
functions (i.e., functions / such that /(xi, . . . , a;„) = 1 — /(I — xi, . . . , 1 — x„)). 

• Graph minor classes (see, e.g., [H]): We assume that the graph class Q is closed 
under taking minors, i.e., G is closed under vertex deletions, edge deletions, and edge 
contractions. Equivalently, Q can be characterized by a finite set of forbidden minors 
[26j . An example is the class of planar graphs (together with the forbidden minors 
i^3,3 and K^). 
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Certainly, other schemes can be devised for system classifications. In fact, many results in 
literature do not fit into this framework (see the discussion of related work below). However, 
the proposed scheme has strong features: first, it exhausts the class of all boolean dynamical 
systems; second, for fixed function classes T and fixed graph classes it is decidable in 
polynomial time whether a given system is in fact, an (.F, ^)-system (supposed the local 
transition functions are represented by lookup-tables); and third, it allows elegant proofs 
of dichotomy theorems (as exemplified by our second contribution). We mention that, as 
fixed points are invariant under changes of the update regime, a scheme for classifying 
update schedules is not needed for our study. 

The main contribution of the paper is a complete complexity classification of the fixed- 
point existence problem with respect to our analysis framework. We make a distinction of 
the problem into three categories: systems for which the local transition functions are given 
by lookup-tables, by formulas (over logical bases), or by circuits (over logical bases). For 
each case we obtain a tractability /intractability dichotomy. Interestingly, the dichotomy 
theorems for formulas and circuits coincide. Let ^ be a Post class of boolean functions 
and let Q be graph class closed under taking minors. The following is proved in Sect. [U 

• Dichotomy for boolean dynamical systems based on lookup-tables: If contains the 
self-dual functions and Q contains the planar graphs, then the fixed-point existence 
problem for (^, ^)-systems with local transition functions given by lookup-tables is 
NP-complete; in all other cases, it is decidable in polynomial time (Theorem 14. ip . 

• Dichotomy for boolean dynamical systems based on formulas or circuits: If J- contains 
the self-dual functions and G contains the graphs having vertex covers of size one, then 
the fixed-point existence problem for (J^, ^)-systems with local transition functions 
given by formulas/circuits (over the logical basis of J-) is NP-complete; in all other 
cases, it is decidable in polynomial time (Theorem 14. Sp . 

The results provide easy criteria for deciding whether (J^, ^)-systems have tractable or 
intractable fixed-point existence problems. For instance, tractability follows for systems 
with linear or monotone local transition function on arbitrary networks (see also [8]). The 
tractability regions with respect to arbitrary local transition functions and restricted graphs 
correspond to bounded treewidth (in the case of lookup-tables) and bounded degree (in the 
case of formulas or circuits). Tractable network classes (for lookup-table-based systems) 
are, e.g., trees, outer-planar graphs, or series-parallel graphs. 

Related work. There is a series of work regarding the complexity of certain computational 
problems for discrete dynamical systems (see, e.g., [171 ETJ El El EJ EQl |28l E] and the 
references therein). 

Detailed studies of computational problems related to fixed-point existence have been 
reported in [8l [30l [28l I29j . As shown in [8j, tractable cases for fixed-point existence are 
constituted by systems with linear, generalized local transition functions, systems with 
monotone local transition functions and systems where local transition function are com- 
puted by gates (of unbounded fan-in) from {AND, OR, NAND, NOR}; intractable cases 
are boolean dynamical systems having local transition functions computed by gates (of 
unbounded fan-in) from the sets {NAND, XNOR}, {NAND, XOR}, {NOR, XNOR}, or 
{NOR, XOR}. Moreover, in [30 1 128 1 129]. the problem of enumerating fixed points of boolean 
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dynamical systems has been studied: counting the number of fixed points is in general #P- 
complete, even counting the number of fixed points for boolean dynamical systems with 
monotone local transition functions over planar bipartite graphs or over uniformly sparse 
graphs is #P-complete. We note that all system classes considered are based on formula 
or circuit representations, i.e., the intractability results fall into the scope of Theorem 14.81 

Mainly, tractability and intractability results have been shown for various version of 
pattern reachability problems such as garden of Eden existence (e.g., [U [30l [28]), pre- 
decessor existence (e.g., [271 [S]), parameterized and unparameterized reachability (e.g., 
[iZl [27l [HJ [U [7] ) . To summarize, the system subclasses considered restrict the local transi- 
tion functions to linear functions, monotone functions, various types of threshold functions, 
or symmetric functions. Except the linear and monotone functions none of these classes is 
closed under composition. Restrictions to the dependency networks involve planar graphs, 
regular graphs, bounded-degree graphs, star networks, and bounded pathwidth. 

In the theory of finite discrete dynamical systems, tight dichotomy results have been 
shown for very restricted classes of systems (as, e.g., in [12]). Exhausting results similar 
to those of this paper are standard for constraint satisfaction problems (see, e.g., [H] for a 
survey). In fact, our results rely on close relationships between fixed-point and constraint 
satisfaction problems (as in the proof of Theorem 14. 3p . 

2 The Dynamical Systems Framework 

In this section, we describe our model of dynamical systems. We follow the approach given 
by [3 with a marginal generalization regarding update schedules. 

Dynamical systems. The underlying network structure of a dynamical system is given by 
an undirected graph G = {V, E) without multi-edges and loops. We suppose that the set 
V of vertices is ordered. So, without loss of generality, we assume F = {1, 2, . . . , n}. For 
any vertex set U , let Ng{U) denote the neighbors of U in G, i.e., 

Ng{U) =dcf { j \ j ^ U and there is an z e [/ such that {i,j} G E }. 

li U = {i} for some vertex i, then we use Ndi) as a shorthand for Nciii})- Define 
Nq^i) =def {i} U Ncii)- The degree di of a vertex i is the number of its neighbors, i.e., 
=dcf \\NG{i)\\. 

A dynamical system S over a domain D is a pair (G, F) where G = {V, E) is an 
undirected graph (the network) and F = {fi \ i E V} is a set of local transition functions 
fi : P'^^+i — >• P. The intuition of the definition is that each vertex i corresponds to an active 
element (entity, agent, actor etc.) which is always in some state Xi and which is capable to 
change its state, if necessary. The domain of S formalizes the set of possible states of all 
vertices of the network, i.e., for all i £ V, it always holds that Xj e P. A vector x = (xj)jgy 
such that Xi € T) for all i € ^ is called a configuration of S. If it is more convenient, then 
we also say that a mapping I : V ^ V is a configuration. A subconfiguration of / with 
respect to j4 C y is a mapping I[A] : A ^ D such that I[A](z) = I{i) for all i £ A. The 
local transition function fi for some vertex i describes how i changes its state depending 
on the states of its neighbors Ncii) in the network and its own state. 

Discrete dynamical systems. We are particularly interested in dynamical system operating 
on a discrete time-scale. A discrete dynamical system S = {S, a) consists of a dynamical 
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system S and a mapping a : {1, . . . , T} — > ViV), where F is a set of vertices of the network 
of S and T € IN. The mapping a is called the update schedule and specifies which state 
updates are realized at certain time-steps: for f € {1, . . . ,T}, a{t) specifies those vertices 
that simultaneously update their states in step t. 

Global maps. A discrete dynamical system S = {S, a) over domain V induces a global map 
Fs '■ T)"" — > where n is the number of vertices of S. For each vertex i £ V, define an 
activity function cpi for a set C/ C.V and x = (xi, . . . , x„) G V^'' by 



where {ii,i2, ■ ■ ■ ,idi+i} = -^g'(^)- ^ U ^ V, define the global transition function 
¥s[U] : ^ for ah x^Vhy 



Note that the global transition function does not refer to the update schedule, i.e., it only 
depends on the dynamical system 5 and not on S. The function Fs : 2^" — > 2?" computed 
by the discrete dynamical system 5, the global map of S, is defined by 



Fixed points. The central notion for our study of dynamical systems is the concept of a 
fixed point, i.e., a configuration which does not change under any global behavior of the 
system. Let S = (G, {fi \ i G V}) be a dynamical system over domain V. A configuration 
X G is said to be a local fixed point of S for [/ C F if and only if Ys\U]{x) = x. A 
configuration x G 2?" is said to be a fixed point of S if and only if x is a local fixed point 
of S for V. 

The following useful proposition is easily seen. 
Proposition 2.1. Let S = {G, {fi \ i G V}) be a dynamical system over domain V. 

1. If the configuration x G is a local fixed point of S for U' ^ V and x is a local 
fixed point of S for U" C V, then x is a local fixed point of S for U' U U" . 

2. A configuration x G is a local fixed point of S for U if and only if x is a local 
fixed point for all U' (^U. 

Notice that the second item of Proposition 12.11 shows that the concept of fixed points is 
independent of update schedules. 

Corollary 2.2. Let S = (G, {fi \ i G V}) he a dynamical system over domain D. A 
configuration x G is a fixed point of S if and only if for all update schedules a : 
{1, . . . ,r} ^ V(y), it holds that Fi^s,a){x) = x. 




Fs[U]{x) =def {ipi[U]{x),...,ipn[U]{x)). 



T 



Fs =deillFs[aik)]. 



k=l 
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3 The Analysis Framework 



In this section, we give a formal description of classification schemes for dynamical systems. 
Local transition functions are classified by Post classes (i.e., closures under composition) 
and networks are classified by forbidden minors (i.e., closures under taking minors). 

Post classes. We adopt notation from [llj. An n-ary boolean function / is a mapping 
/ : {0, 1}" {0, 1}. Let BF denote the class of all boolean functions. There are two 1-ary 
boolean functions: id(3;) =def x and not(x) =dcf 1 — x (which are denoted in formulas by x 
for id(x) and x for not(x)). For b G {0, 1}, a boolean function / is said to be b-reproducing 
if and only if f{b, . . . ,b) = b. For binary n-tuples a = (oi, . . . , o„) and b = (bi, . . . , bn), we 
say that (oi, . . . , a„) < (bi, . . . , 5„) if and only if for all i € {1, . . . , n}, it holds that < bi. 
An n-ary boolean function / is said to be monotone if and only if for all x, y G {0, 1}"', 
X < y implies f{x) < f{y). An n-ary boolean function / is said to be self-dual if and only 
if for all (xi, . . . , Xn) € {0, 1}", it holds that f{xi, . . . , x„) = not(/(not(xi), . . . , not{xn))). 
A boolean function / is linear if and only if there are constants ai, . . . , a„ G {0, 1} such 
that f{xi, . . . , Xn) = ao © oixi © • • • © a„x„. Here, © is understood as addition modulo 2 
and xy is understood as multiplication modulo 2. 

We say that a function class J- is Post if and only if contains the function id and 
is closed under introduction of fictive variables, permutations of variables, identification of 
variables, and superposition (see, e.g., [11] for a formal definition). It is a famous theorem 
by Post [23] that the family of all Post classes of boolean functions is a countable lattice 
with respect to set inclusion (see, e.g., [231 |3l] a proof). The maximal meet-irreducible 
classes are the following (see, e.g., [TT]): 



Ro 


=def 


{ / G BF 


/ 


is 


0-reproducing } 


with logical basis {AND,XOR} 


Ri 


=def 


{ / G BF 


/ 


is 


1-reproducing } 


with logical basis {OR, XNOR} 


L 


=def 


{ / G BF 


/ 


is 


linear } 


with logical basis {XOR, 0, 1} 


M 


=def 


{ / G BF 


/ 


is 


monotone } 


with logical basis {AND, OR, 0, 1} 


D 


=dcf 


{ / G BF 


/ 


is 


self-dual } 


with logical basis 



{{x Ay) V (x A z) V {y Az)} 



From the structure of Post's lattice follows that each other class which does not contain 
all boolean functions is included in an intersection of two of the classes listed. 

Graph minor classes. We adopt notation from [TJ]. Let X and Y be two undirected graphs. 
We say that X is minor of Y if and only if there is an isomorphic subgraph Y' of Y such that 
X is obtained by contracting edges of Y' . Let ^ be the relation on graphs defined hj X ^Y 
if and only if X is a minor of Y. A class G of graphs is said to be closed under taking minors 
if and only if for all graphs G and G' , if G G ^ and G' ^ G, then G' G G- Let X be any set of 
graphs. Forb^ {X) denotes the class of all graphs without a minor in X (and which is closed 
under isomorphisms). More specifically, Forb-<(^) =dcf {G \ G ^ X for all X G ^ }. The 
set X is called the set of forbidden minors. Note that Forb^(0) is the class of all graphs. 
As usual, we write Forb^(Xi, . . . , Xn) instead of Forb^({Xi, . . . , Xn}). A useful property 
of the forbidden-minor classes is the monotonicity with respect to ^, i.e., X ^Y implies 
Forb^(A:) C Forb^(y). 

The celebrated Graph Minor Theorem, due to Robertson and Seymour [26], shows 
that there are only countably many graph classes closed under taking minors: A class G 
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of graphs is closed under taking minors if and only if there is a finite set X such that 
Q = Forb^(Af). An important consequence of this theorem is that for each graph class 
Q which is closed under taking minors there exists a polynomial-time algorithm to decide 
whether a given graph belongs to G [25]. 

The most prominent example of a characterization of a class closed under taking minors 
in terms of forbidden sets are the planar graphs. Let K"^ denote the complete graphs on n 
vertices and let Kn^m denote the complete bipartite graph having n vertices in one compo- 
nent and m vertices in the other component. The well-known Kuratowski- Wagner theorem 
(see, e.g., [M]) states that a graph G is planar if and only if G belongs to Forb^(-ftr3^3, K^). 
Planar graphs have an algorithmically important property: a graph X is planar if and 
only if Forb-<(X) has bounded treewidth [23|. As we use treewidth of a graph only in a 
black-box fashion, we refer to, e.g., [H] for a definiton of treewidth. A class Q of graphs is 
said to have bounded treewidth if and only if there is a A; G IN such that all graphs in the 
class have treewidth at most k. 

A similar but much less subtle behavior to planar graphs show graphs with a vertex 
cover of size one. Let G = {V, E) be a graph. We say that a subset {/ C 1/ is a vertex cover 
of G if and only if for all edges {u, f } G it holds that {n, n C/ ^0. It is known that the 
class of graphs having a vertex cover of size at most k is closed under taking minors |13j . 
Moreover, G has a vertex cover of size one if and only if G belongs to Forb-< (X^, © K"^) 
|13j . where for graphs G and G', G (B G' denotes the graph obtained by the disjoint union 
of G and G' . A class of graphs is said to have bounded degree if and only if there is a A; € IN 
such that all graphs in the class have a maximum vertex-degree of at most k. 

Proposition 3.1. Let X be a graph. Then, X has a vertex cover of size one if and only 
i/Forb^(X) has bounded degree. 

Proof. For (^), suppose that X has a vertex cover of size one. Then, X consists of some 
star graph Ki^k and some isolated vertices ui,...,Ur. Assume that Forb^(X) does not 
have bounded degree. Thus, there exists a graph G € Forb-<(X) of maximum vertex- 
degree at least k + r. So, G contains a subgraph Ki^k+r- Hence, X ^ G, contradicting 
G € Forb^(X). Therefore, Forb^(X) has bounded degree. 

For (<^=), suppose that X does not have any vertex cover of size one. It is easily 
seen that in this case, X contains a triangle or two non-incident edges. First, suppose X 
contains a triangle K^. As Foih^{K^) contains the class of all trees, which certainly does 
not have bounded degree, we easily obtain from ^ X, that Forb-<(X) does not have 
bounded degree as well. Second, suppose X contains at least two non- incident edges, i.e., 
K'^®K'^ < X. It follows that Forb^(X) contains for ah /c G IN, the star graph Ki^^. Thus, 
Forb-< (X) does not have bounded degree. This completes the proof of the direction from 
left to right. □ 

4 Islands of Tractability for Fixed-Point Existence 

In this section we are interested in the computational complexity of the following problem. 
Let .F be a class of boolean functions and let ^ be a class of graphs. 
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Problem: 
Input: 



Question 



FixedPoints(:F, g) 

A boolean dynamical system S = (G, {/i, . . . , /„}) such that G G G and 
for alH G {1, . . . , n}, fi^T 
Does S have a fixed point? 



The complexity of the problem depends on how transition functions are represented. We 
consider the cases of look-up table, formula, and circuit representations. The correspond- 
ing problems are denoted by FixedPointSt, FixedPointSf, and FixedPointSc. It is 
obvious that all problem versions belong to NP. We say that a problem is intractable if it 
is NP-hard, and it is tractable if it is solvable in polynomial time. 

4.1 The Case of Local Transition Functions Given by Look-up Tables 

The main result of this subsection is the following dichotomy result. 

Theorem 4.1. Let T he a Post class of boolean functions and let Q he a class of graphs 
closed under taking minors. IfT andg D Forh^iKs^s, K^) i/ien FixedPointSt(J^, ^) 
is intractable. Otherwise, FixedPointSt(.^, ^) is tractable. 

We postpone the proof to the end of this paragraph when we have proved a number of 
complexity results for several classes of systems. 

Dynamical systems with a tractable fixed-point analysis. We first summarize the tractable 
classes of local transition functions of systems on arbitrary underlying networks. The 
results presented next are well-known or follow easily from the definitions. 

Proposition 4.2. 1. FixedPointSt(Ri, Forb(0)) is solvable in polynomial time. 

2. FixedPointSt(Roi Forb(0)) is solvable in polynomial time. 

3. [8] FixedPointSt(M, Forb(0)) is solvable in polynomial time. 

4- [8] FixedPointSt(L, Forb(0)) is solvable in polynomial time. 

Notice that D is the only remaining class which does not contain all boolean functions. 

Next we restrict network classes. To identify tractable cases, we express the fixed- 
point existence problem as a constraint satisfaction problem. A constraint satisfaction 
problem (CSP) consists of triples {X, V, C), where X = {xi, . . . , x„} is the set of variables, 
T> is the domain of the variables, C is a set of constraints Rxi^, . . . ,Xi^ having associated 
corresponding relations Ri^^...^if., such that there exists an assignment I : X ^ V satisfying 
),..., € Rij^,...,i^^ for all constraints iix^^ , . . . , Xi^. € C. We suppose that C is 
listed by pairs {Rxi^,. . . , Xj^, iiij,...,^^^). 

Theorem 4.3. Let X be a planar graph. Then, FixedPointSt(BF, Forb^(X)) is solvable 
in polynomial time. 

Proof. We describe a general reduction from FixedPointSt(.^, ^) to constraint satisfac- 
tion problems. Suppose we are given a dynamical system S = (G, {/i, . . . , /„}). Let 
G = {V^ E) be the underlying network. Without loss of generality, we may assume that G 
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does not have isolated vertices, i.e., dj > 1 for all i &V. Define CSP(S') = {X,V,C) to be 
the constraint satisfaction problem specified as follows: 

X =def {xi,...,Xn} 

=def UieV '^i where for all i & V, 

Vi =def {(/, ^) I / : N^ai^ ^ {0, 1} such that /^(/(io), • • • , I{ik)) = Hi) 
where {io, ...,ik} = N^{i)} 

C =def {ExiXj I {i,j} G E{G) and i < j} where for all i < j, 

E^J =def { {{lui), I Ii{k) = Ij{k) for all k G N^{i) n N^{j) } 

Let n = and m = By construction of CSP(S') we obtain that ||X|| = n and that 
the number of constrains in C is just m. The size of the domain P is at most proportional 
to EievC^ + di) ■ 2^+'^% and the size of the set of constraint relations can bounded by 

for some constant c > 0. The latter holds because di > 1 for all i G V. All in all, this 
easily implies that |CSP(S')[ = 0{\S\'^). Hence, CSP(5) is computable in polynomial time 
in the size of S. We have to show that 

S has a fixed point <^=^ CSP(S') has a satisfying assignment. 

We prove both directions of the equivalence separately. 

For (=^), suppose the configuration I : V ^ {0,1} is a fixed point of S. Define an 
assigment /' : {xi, . . . ,Xn} — * by I\xi) =def 

I[NQ{i)\. We show that /' satisfies all 
constraints. For all i £ V, let Ij denote I[NQ[i)]. Let ExiXj be any constraint of CSP(S') 
and let Eij be the relation associated with ExiXj. By definition, € E. Since / is 

a fixed point, it holds that and belong to V. For all k E ^^{i) n NQ{j) we 

obtain Ii{k) = I'{xi){k) = I'{xj){k) = Ij{k). Thus, lies in Eij. This proves 

the direction from left to right. 

For (^), suppose / is a satisfying assignment for CSP(S'). Let £ V he the 

pair which / assigns to the variable Xi G V. Define a configuration I' : V ^ {0, 1} by 
I'ii) =dcf liii)- We are done if we are able to show that the following is true for all i (zV: 

l'[NUi)]^h (4.1) 

Since for all i £ V, fi{Ii{iQ), . . . , Ii{ik)) = where {io,...,'jfc} = Eq. (fiTll 

implies that /' is a fixed point of S. To verify the equation, we consider i and all neighbors 
of i individually. Let j E NQ^i). If j = i, there is nothing to show. So let j ^ i- 
Since G E{G) and since / is a satisfying assignment, Ii{k) = Ij{k) holds for all 

k € -^G'(^) ^G^j)- Thus, we obtain = = hij)- This shows the correctness of 
Eq. ()4.ip for all i €V. Hence, the direction from right to left is proved. 

We thus have established a polynomial-time many-one reduction between the fixed- 
point existence problem and constraint satisfaction problems. Moreover, it is easy to see 
that the graph G of a given dynamical system S is isomorphic to CSP(S')'s constraint 
graph which consists of the vertex set {xi, . . . , Xn} and the edge set {{xi,Xj} \ ExiXj € C}. 
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Thus, if X is planar, then Forb^ (X) has bounded treewidth, and so the constraint graph of 
CSP(5') has bounded treewidth for each dynamical system S = {G, {/i, . . . , /„}) such that 
G G Forb-<(X). The theorem follows from the well-known polynomial-time algorithms for 
constraint satisfaction problems having constraint graphs of bounded treewidth [15]. □ 

Dynamical systems with an intractable fixed-point analysis. We turn to the intractable 
cases of the fixed-point existence problem. As a first step, FixedPointSt(.^, ^) is shown 
to be NP-complete for arbitrary boolean, local transition functions and planar networks. 
Recall that similar results have only been shown for formula respresentation. 

We use Planar 3SAT to prove NP-hardness. Let H = Ci A - • - A Cm be a propositional 
formula in conjuctive normal form where each clause Cj consists of three literals (for short, 
3CNF) where positive literals Xi and negative literals -iXj are taken from the set of variables 
{xi, . . . ,Xn}- The graph representation T{H) of -ff is a bipartite graph consisting of the 
vertex set {xi, . . . ,Xn,Ci, . . . , Cm} and all edges {xi, Cj} such that variable Xi appears as 
a literal in the clause Cj. A planar 3CNF is a 3CNF such that its graph representation is 
planar. Planar 3SAT is the problem to decide whether a given planar 3CNF is satisfiable. 
It is well known that Planar 3SAT is an NP-complete problem |19j . 

Theorem 4.4. FixedPointSt(BF, Forb(K3^3, K^)) is ^Y" -complete, even restricted to un- 
derlying networks having maximum vertex degree three. 

Note that for graphs having maximum degree at most two the problem is tractable. 

Proof. Let if = Ci AC2A - • - A Cm be a 3CNF having variables Xi , . . . , Xn and a planar graph 
representation T{H) where V{T{H)) = U [JW , U = {xi, . . . and W = {Ci,.. . ,Cm}. 

Suppose that V{r{H)) is totally ordered such that the clause vertices of W come before 
the variable vertices of U . 

We construct the following system. For the underlying network C = (y,E), compute 
an embedding of T{H) in the plane (in linear time). Now replace the vertices of U in 
the following way. Let Xi G U and suppose that Xj's neighbors Cj-^, . . . ,Cj^. in T{H) 
are clockwise ordered with respect to the planar embedding, such that j'l is minimal in 
{j I Xi appears as a literal in Cj }. Replace the vertex x, € C/ by a cycle {xjj^, . . . ,Xij^} 
which is connected to the neighbors of Xj in T{H) by having an edge {xij, Cj} whenever Xj 
appears as a literal in Cj. Let C = {V, E) be the graph obtained after all replacements are 
made in this way. Clearly, G is planar, i.e., G G Forb^(-K'3^3, ii'^), and can be computed 
in time polynomial in the size of H. Note that the maximum degree of a vertex in C is 
three. To complete the construction, we specify the local transition functions. Let Ci gW 
be clause vertex. Let a^n ,ji , a^jj j2 ) ^-^d Xi^^j^ be the neighbors of Cj in C. Define the local 
transition function /c- by 

fc, (liCi), I{Xi^ ) , I(Xi2 ,32 ) . li^h ,j3 ) ) 

1 if {l{xi^j^), I{xi2j2)j H^isJs)) is a satisfying assignment of C, 

not(I(Cj)) otherwise 



-def 



Let Xij G y be a variable vertex. Suppose that {Cj, Xj^fcoj • • • ) = ^ci^ij) where r is 

the degree of Xij in G. Define the local transition function f^^ ^ by 



-dof 



I{xij) if I(xi,fco) = • • • = HHkr) 

not(/(xjj)) otherwise 
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Let Sh denote the dynamical system (G, € V} constructed from any planar 3CNF 
H in the way just specified. Note that, since the maximum degree of any vertex in G is 
three, we can compute Sh in time polynomial in the size of H. 

By construction of Sh, a configuration / is a fixed point of Sh if and only if I satisfies 
H^iji) = ••• = I{xij^) for all i € {!,..., n} and, furthermore, I{Cj) = 1 for all j G 
{!,... ,m}. Hence, it is easily seen that H has a satisfying assignment if and only if Sh 
has a fixed point. □ 

We extend the NP-completeness of the fixed-point existence problem to systems with 
self-dual transition functions and planar graphs. First, observe that the transition func- 
tions for clause variables are not self-dual. This implies that the construction used in the 
last theorem is not fully appropriate. However, arbitrary boolean functions can easily be 
embedded into a self-dual function of larger arity. 

Proposition 4.5. Let n G IN+. For each k-ary boolean function f : {0, l}'^ {0, 1}, the 
function sd„(/) : {0, Ij'^+^+i {0, 1} defined for all xi, . . . , Xk, yi, ■ ■ ■ ,yn,z e {0, 1} by 

{fixi,...,Xk) ify^ = ... = y^ = 

not(/(not(2;i), . . . ,not(xfc))) if y^ = ■ ■ ■ = y^ = I 
not(z) otherwise 

is self-dual. 

Proof. Follows from the definitions by case analysis. □ 

The usage of Proposition 14.51 introduces ambiguity to the set of fixed points. 

Proposition 4.6. Let S = {G, {fi \ i G V}) be a dynamical system over {0, 1} so that all 
local transition functions fi are self-dual. Let [/ C y. Then, a configuration / : 1/ — > {0, 1} 
is a local fixed point of S for U if and only if the configuration / : 1/ — > {0, 1} defined by 
I{i) = not(/(i)), is a local fixed point of S for U. 

Proof. Immediate from the definitions. □ 
Theorem 4.7. FixedPointSt(D, Forb(iir3^3, K^)) is NF- complete. 

Proof We reduce from FixedPointSt(BF, Forb(i^3,3, i^^)). Let S = {G,{fi \ i G V}) be 
a dynamical system such that the underlying network G is planar. We construct another 
system S' = | i G V'}) as follows. The network G' = {V',E') is defined by 

V =def V U E and E' =def E U { {i,e} \ i £ V,e € E,i e e }. It is easily seen that G' is 
planar as well. Suppose V' is ordered such that V is in the same ordering as for S, E is 
arbitrarily ordered, and E comes completely after V. The local transition functions of the 
vertices of V are specified as follows. Suppose i € V. Let {io, . . . ,ik} = -^c-'(^) ^ 
let {ei, . . . , Cfc} = N^,{i) n E. We define the local transition function f^ by 

fi {"^io 1 • • • ) "^ik ' '^^1 ' • • • ) ■^Gfe ) — def sdfc (/j) (jjjg , . . . , ^^j^. , Xgj , . . . , 3^6^ j ) • 

By Proposition 14.51 and since D is closed under identification of variables, /| is self-dual. 
Moreover, as the degree of i in V' is 2k where k is the degree of i in G, the function table 
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can be computed in polynomial time. Now consider a vertex e = {i,j} G E. Define the 
local transition function /g by 

Clearly, is self-dual. As the degree of e € -E is two, the function table is trivially 
computable in polynomial time. It remains to show that 

5 has a fixed point <;=^ S' has a fixed point. 

We prove both direction individually. 

For (=^), suppose the configuration / : 1/ — > {0, 1} is a fixed point of 5, i.e., for all i € V, 
it holds that /j(/(io), • • • , I{ik)) = where {io, . . . ,ik} = ^^{i). Define a configuration 
I' -.V ^ {0, 1} for iGVhj =dcf and for e G S by I'(e) =def 0. Consider i e V. 
Assume {io, • • • , 4} = ^g"(^) ^ ^""^ i^i' • • • ' '^k] = ^g'(^) ^ ^- obtain 
/;(/'(io),...,/'(ife),/'(ei),...,/'(efc)) 
= sdfc(/i)(/(io),...,/(4),0,...,0,/(f)) = /i(/(ii),...,I(4)) = /(i) = /'«. 

Thus, /' is a local fixed point for i ^ V. Suppose e = G -B. By definition, 

f[^ .^{I'{i)J'{j)J'{{i,j})) = r{{i,j]). Thus, /' is a local fixed point for e^E. Proposi- 
tion [2T] implies that /' is a fixed point of S' . 

For (<^=), suppose the configuration I' : — > {0, 1} is a fixed point of S' . Observe that 
for all e, e' G C y , there is a walk in G' from e to e' alternating between vertices in V 
and i?, i.e., there are vertices po-,. ■ ■ ,P2e G V such that {pi,pi^i} G i?' for all < z < 2i, 
Po = e, = e', and for all < j < it holds that G E and P2j+i G Consider a 
vertex p2i+i G V. Let {zq, • • • , ifc} = N^,{p2j+i) n 1/ and {ei, . . . , Cfc} = N^,{p2j+i) n 
Since /' is a fixed point of S, we have 

/'(P2,+i) = 4^.^,(/'(io),...,/'(4),/'(ei),...,/'(efc)) 

= sdfc(/p,^.,, (/'(io), . . . , /'(ifc), I'(ei), . . . , /'(efc), /'(i)) 

By definition of sdk{fp2j+i), it follows that I'{ei) = • • • = /'(cfc). In particular, I'{p2j) = 
I'{p2j+2)- This implies that for all e,e' G -E, it holds that /'(e) = I'{e'). By Proposition 
14.61 we may assume that /'(e) = for all e £ E. Define a configuration I : V ^ {0, 1} of 
S by /(i) = /(i) for all i £V. It is easily seen that / is a fixed point of S. □ 

Composing the big picture. We come back to the proof of the main result of the subsection. 
For convenience, we state it once more. 

Theorem 14.11 Let T he a Post class of boolean functions and let Q be a class of graphs 
closed under taking minors. IfJ^ andg D Forb-<(/C3,3, /C^) t/ien FixedPointSt(J^, ^) 
is intractable. Otherwise, FixedPointSt(.^, ^) is tractable. 

Proof If 2 D and ^ 2 Forb(/r3,3, /C^), then FixedPoints(J', ^) is NP-complete by 
Theorem 14.71 Suppose J- ^ D or Q ^ Forb(/C3^3, /C^). The maximal classes J- that do 
not contain D are Ri, Rq, M, and L. For all these classes, by Proposition 14.21 the fixed- 
point existence problem is solvable in polynomial time. It remains to consider the case 
that g 2 Forb(/C3,3,i^^)- Suppose G = Forb(Xi, . . . , X^). Since G 2 Forb^(/C3,3, ^^), 
there is an i such that Xi is planar. Since G ^ Forb-<(Xj), Theorem 14.31 shows that 
FixedPointSt(BF, G) is solvable in polynomial time. □ 
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4.2 Succinctly Represented Local Transition Functions 

In this section we prove a dichotomy theorem for the fixed-point existence problem when 
transitions are given by formulas or circuits. As usual, the size of formula is the number of 
symbols from the basis used to encode the formula, the size of a circuit is the number of 
gates it consists of (including the input gates). Both succinct representations of functions 
lead to the same result. 

Theorem 4.8. Let T he a Post class of boolean functions and let Q be a class of graphs 
closed under taking minors. IfJ^ andg ^ Forh^{K^ , K"^ ^ K^) , then FixedPointSf 
{J^,Q) is intractable. Otherwise, FixedPointSf(.?^, ^) is tractable. Moreover, the same 
statement holds for FixedPointSc. 

Again we postpone the proof until we have proved a number of special results. 

On the side of tractable cases, first note that Proposition 14.21 still holds for formulas 
and circuits. Actually, the results in [8] were stated for formulas. Furthermore, notice 
that circuits over the basis {XOR, 0, 1} can easily be transformed in polynomial time into 
equivalent formulas over the same basis. The following result provides the tractability limit 
for restricted network classes. 

Theorem 4.9. Let X be a graph having a vertex cover of size one. Then, FixedPointSf 
(BF, Forb^(X)) and FixedPointSc(BF, Forb^(X)) are solvable in polynomial time. 

Proof. Suppose X has a vertex cover of size one. Then, there is an r S IN such that all 
graphs in Forb-<(X) have maximum vertex degree r. Thus, if we compute lookup tables 
from formulas or circuits, then each lookup table has at most 2^~^^ entries. Hence, in 
polynomial time, we can transform each dynamical system 5 with a network in Forb-<(X) 
and local transitions functions given by formulas or circuits into a dynamical system S' 
with the same networks and local transition functions given by lookup tables such that 
S and S' have the same fixed-point configurations. Moreover, since X is planar (note 
that ^ ^3,3 and K'^ ® K"^ ^ Ks^3 as well as ^ i^3,3 and K'^ ® K'^ < K^), 
Theorem 14.31 imphes polynomial-time solvability of FixedPointSf(BF, Forb^(X)) and 
FlXEDPoiNTSc(BF,Forb^(X)). □ 

Theorem 4.10. FixedPointSf(D, Forb^(i^^, © ET^)) is NF- complete. 

Proof. The reduction is from 3SAT. We start with a description of a reduction to dynami- 
cal systems where each transition function is computed by a 3CNF. Suppose we are given a 
3CNF H = CiA- ■ -ACm having variables xi, . . . ,Xn. Note that the formula H' =(iof -ff Vxq, 
where xq is a new variable, satisfies for any assignment / : {xq, xi, . . . , x„} {0,1} 
that fH'{I{xo), I{xi), . . . ,L{xn)) = I{xo) if and only if I{H) = 1 and I{xq) = 1. More- 
over, H' is equivalent to a 4CNF which can be transformed in a 3CNF H with vari- 
ables xq,xi, . . . ,Xn, Xn+i., • • • , x„_|_m such that H is satisfiable with 1 assigned to xq if and 
only if H is satisfiable with 1 assigned to xq. Define Sh to be the dynamical system 
{G, {/o, . . . , /n+m}) consisting of the network G = {V, E) where V =Aei {0, 1, . . . , n + m} 
and E =Aci { {0,i} [ i € + m} }. Thus, G is a star Ki^n+m-, i-e., G € 

Forb^(i^^, ii'^ K"^). The local transition functions are given as follows. For a vertex 
i € {1, . . . , n + m}, the local transition function /j is defined to be computed by the for- 
mula i?(xo,Xj) =def Xi. For the central vertex, the local transition function /q is given by 
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the formula H{xo,xi, . . . ,Xn+m) where the variable Xi stands for a vertex i G V. Clearly, 
Sh can be computed in time polynomial in the size of H and we have that H is satisfiable 
if and only if Sh has a fixed point. 

We now transform the dynamical system Sh into another system S'jj with self-dual 
local transition functions given by formulas over the corresponding basis. Note that D has 
a single basis function of arity three. Let D denote the corresponding ternary function 
symbol, i.e., the semantics of D is defined by D(x, y, z) =dof (x Ay) V {x A'z) V {y A'z). Note 
that D(x,x,y) = y. We embed a 3CNF into a self-dual function, similarly to Proposition 
14. 5i That is, for an arbitrary 3CNF H = Ci A • • • A Cm having variables xi,...,Xn, 
define the formula dual(i7)(xi, . . . , x„) = Hixi, . . . , x^). Define sd(//)(xi, . . . , x„, z) =def 
{H A V (dual(i?) A l). By induction over the number of clauses, we show that s(i{H) is 
equivalent to a formula built using D which is of polynomial size: 

1. For the base of induction, suppose m = 1. Since we know how to express negation 
using L), without loss of generality, we assume that H = {xi V a;2 V xs). So, sd(i/) = 
((xi V V X3) A V ((xi Ax2 A x^) A z) . By truth-table inspection we obtain that 
sd(i?) = D{z,D{z,xl,x^),D{z,xl,x^)). 

2. For the induction step, suppose m > 1. Let H' = Ci A ■ ■ ■ A Cm be a 3CNF over 
the variables xi, . . . , x„. Define H' =Aci Ci A • • • A C\m/2\ ^-nd H" =Aei C[m/2j+i ^ 
■ ■ ■ A Cm- Some equivalent transformations show that sd[H) = D{'z, sd{H'),sd{H")). 
By induction hypothesis, sd(H') and sd{H") can be expressed using D. Replacing 
sd{H') and sd{H") (and the negation) gives the appropriate formula for sd{H). Note 
that the recursion depth for formula replacement is logarithmic. It follows that the 
size of the formula for sd(-ff) is 0{\H\'^). 

Finally, we define the dynamical system S'jj to be specified as follows. The network C = 
{V, E') consists of the vertex set V =dei {0, 0, . . . , n + m} and the edge set E' =Aei 
{{0,i} I i G {0, ... ,n + m} }. Thus, C is a star i.e., C G Forb-<(K^ e K^). 

The local transition function /j for a vertex i G {0, . . . ,n + m} is given by H{xQ,Xi) = 
D{xo, xq, D{xo, xq, Xi)) . Notice that H{xQ,Xi) = Xj. The local transition function /o for 
the vertex is given as follows. Recall that H is the 3CNF associated with the local 
transition function of vertex in the system Sh ■ Then /o is represented by the L)-formula 
equivalent to sd{H){xo, . . . ,Xn+m,xo). Clearly, S'^ can be computed in time polynomial 
in the size of H. Moreover, it is easy to verify that H is satisfiable if and only if S"^ has a 
fixed point. □ 

We combine Theorem 14.91 and Theorem 14. 101 to prove Theorem 14.81 

Proof. (Theorem\4^ If D D and ^ 2 Fovh{K^,K'^ K"^), then FlXEDPoiNTS(:r, ^) 
is NP-complete by Theorem HjOl Suppose JF ^ D or ^ Forb(K^K2 e K"^). The 
maximal classes J- that do not contain D are Ri, Rq, M, and L. For all these classes, by 
Proposition l4.2l the fixed-point existence problem is solvable in polynomial time. It remains 
to consider the case that G 2 FoTh{K^ , K'^ © K"^). Suppose G = Forb(Xi, . . . Since 
G 2 Forb^(/^^, © K'^), there is an i such that Xi has a vertex cover of size one. Since 
G C Forb^(Xj), Theorem 14.91 shows that FixedPointSt(BF, C/) is solvable in polynomial 
time. □ 
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5 Conclusion 



We characterized the islands of tractabiUty for the fixed-point existence problem for boolean 
dynamical systems with respect to transition classes J- closed under composition and net- 
work classes Q closed under taking minors: \i!F contains the self-dual functions and Q con- 
tains the planar graphs, then FixedPointSt(-?^, ^) is intractable, otherwise it is tractable. 
Replacing "planar graphs" by "graphs having a vertex cover of size one" yields the same 
dichotomy theorem for the succinct representations of local transition functions by formu- 
las or circuits. The linear and monotone functions have been shown to be tractable cases 
in [1]. There, the authors suggested to find more tractable classes of local transition func- 
tions. Over the boolean domain our results show that, aside from two obvious exceptions 
(the 0- and 1-reproducing functions), there are no more such function classes. 

Although the proposed analysis framework allows elegant dichotomy theorems for fixed- 
point existence, it is certainly necessary to examine its usefulness for other computational 
problems for discrete dynamical systems and to refine it appropriately. Another important 
open problem is the extension of the dichotomy theorems to larger domains. This seems 
a challenging issue as even in the case of a ternary domain the number of Post classes is 
continuum (see, e.g., [22]). 

Acknowledgements. I am grateful to Ernst W. Mayr (TU Miinchen) for careful proof- 
reading and pointing out an error in earlier versions of this paper. 
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